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ABSTRACT 



Aims. The nearby galaxy clusters Abell 496 and Abell 85 are studied in the very high energy (VHE, E> 100 GeV) band to 
investigate VHE cosmic rays (CRs) in this class of objects which are the largest gravitationally bound systems in the Universe. 
Methods. H.E.S.S., an array of four Imaging Atmospheric Cherenkov Telescopes (lACT), is used to observe the targets in the 
range of VHE gamma rays. 

Results. No significant gamma-ray signal is found at the respective position of the two clusters with several different source size 
assumptions for each target. In particular, emission regions corresponding to the high density core, to the extension of the entire 
X-ray emission in these clusters, and to the very extended region where the accretion shock is expected, are investigated. Upper 
limits are derived for the gamma-ray flux at energies E> 570 GeV for Abell 496 and E> 460 GeV for Abell 85. 
Conclusions. From the non-detection in VHE gamma rays, upper limits on the total energy of hadronic CRs in the clusters are 
calculated. If the cosmic-ray energy density follows the large scale gas density profile, the limit on the fraction of energy in these 
non-thermal particles with respect to the total thermal energy of the intra-cluster medium (ICM) is 51% for Abell 496 and only 
8% for Abell 85 due to its larger mass and higher gas density. These upper limits are compared with theoretical estimates. They 
predict about ~10% of the thermal energy of the ICM in non-thermal particles. The observations presented here can constrain 
these predictions especially for the case of the Abell 85 cluster. 
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1. Introduction 

Clusters of galaxies, the most extended gravitationally 
bound systems in the Universe, are expected to contain a 
significant population of non-thermal particles. The mech- 
anisms for producing cosmic rays (CR) in galaxy clusters 
can be divided into external processes and internal pro- 
cesses. In external processes, the particle acceleration is 
driven by the assembly of the cluster. An efficient pro- 
duction of high energy particles is expected especially at 
the strong accretion shock at the outskirt of the cluster, 
where cold in-falling material plunges into the already ex- 
isting hot intra-cluster medium (ICM). Therefore, large- 
scale shock waves caused by the cosmological structure 
formation may populate these objects with a non-thermal 
component of particles (see e.g. Colafrancesco & Blasi 
fT998l Loeb & Waxman [20Q01 Ryu et al. [20031 Miniati 
120031) . Particles can also be accelerated by turbulence in 
the intra-cluster medium(ICM) generated by major sub- 
cluster merger events (e.g. Brunetti et al. 12004^ . In con- 
trast in the internal mechanisms, the CRs are accelerated 
by cluster galaxies and injected into the whole cluster vol- 
ume afterwards. Internal sources of CRs in clusters can 
be supernova-driven galactic winds (Volk et al. I1996|) or 
AGNs (e.g. EnBlin et al. lT997l Aharonian l2002l Hinton et 

ai. mm . 

Accelerated, highly energetic CR particles can subse- 
quently produce gamma rays through two main mecha- 
nisms in galaxy clusters. VHE gamma rays can be gener- 
ated by hadronic processes, through inelastic collisions be- 
tween high-energy protons and target protons provided by 
the thermal hot ICM and subsequent 7r°-decay (Dennison 
I1980p . In the leptonic scenario, very high energy (VHE, 
E > 100 GeV) electrons up-scatter cosmic microwave 
background (CMB) photons to gamma-ray energies via 
the inverse Compton process (e.g. Atoyan & Volk 120001 
Gabici & Blasi [20031 [2004]) . VHE electrons which up- 
scatter low energy photons to the gamma-ray range can 
also be of secondary origin either due to photon - photon 
interactions (Timokhin et al. I2004p or generated by ultra- 
high energy protons (E > 10^* eV) interacting with CMB 
photons in the Bethe-Heitler process (Inoue et al. I2005p . 

Galaxy clusters are huge structures with a linear scale 
of several Mpc. This length scale together with typical 
cluster magnetic fields in the /iC range (Carilli & Taylor 
[^00^ Govoni & Feretti [MHj) leads to a large diffusion 
timescale. This diffusion time is longer than the Hubble 
time for CR protons with energies less than '^lO^^ eV. 
Therefore no losses of such particles will occur in these 
systems (Volk et al. 119961 Berezinsky et al. ll997p . Gamma 
rays produced through the hadronic scenario can thus be 
used to probe particle acceleration in clusters of galaxies 
over the entire formation history. 
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In contrast to hadronic CRs, electrons which are accel- 
erated to VHE energies in clusters will suffer from substan- 
tial energy losses due to synchrotron and inverse Compton 
radiation. Therefore, an electron with an energy of ~1 
TeV has only a limited lifetime of ~10^ years in a typical 
cluster magnetic field of 1 fj,G (e.g. Atoyan & Volk 12000]) . 
Consequently VHE electrons do not accumulate in clus- 
ters and it is expected that the electron spectrum will cut 
off at an energy that is below the VHE regime outside the 
acceleration sites. Hence leptonic VHE gamma-ray pro- 
duction is likely less relevant in galaxy clusters (Atoyan & 
Volk l20n0)l . 

Observationally, a significant population of non- 
thermal electrons is found in several clusters that can 
be detected in radio waves (Giovannini & Feretti 120001 
Feretti et al. [20041 Bagchi et al. 2006) and possibly hard X- 
rays (Rephaeh & Gruber|2002l Fusco-Femiano et al [2004|) . 
Given the evidence/hints of accelerated particles in galaxy 
clusters, these objects are potential sources for observation 
of gamma-rays (see Blasi et al. 120071 for a recent review). 
However, despite these indications, no cluster has been 
established as a gamma-ray emitter so far (Reimer et al. 
[2003]) . In the VHE gamma-ray range Perkins et al. (p006|) 
have reported upper limits for two nearby massive clusters 
(Perseus and Abell 2029) with the Whipple telescope. For 
the Perseus cluster, these results imply that the cosmic- 
ray proton energy density is less than 8% of the thermal 
energy density. 

Throughout the paper a standard cosmology with 
r^A = 0.7, rim — 0.3 and Hq = 70 km s~^ Mpc~^ is used. 

In this paper, the observations of the galaxy clusters 
Abell 496 and Abell 85 with the H.E.S.S. array of tele- 
scopes in the VHE gamma-ray regime are presented. 

2. Target overview 

The target clusters were selected in terms of optimal de- 
tectability, position and distance for an observation with 
H.E.S.S. Promising targets of this kind should be located 
on the southern hemisphere and at a redshift not larger 
than z ^ 0.06, since more distant objects suffer substantial 
absorption from extragalactic background light (EBL). 
Furthermore, there should be no blazar at the location 
of the cluster which could superpose potential VHE emis- 
sion of the galaxy cluster. Two of the clusters that were 
chosen for an observation with H.E.S.S. are presented in 
this paper. These clusters were selected according to dif- 
ferent criteria as described below. Theproperties of both 
selected clusters can be seen in Tab. [Jo. 

As a first target, a compact galaxy cluster was selected. 
The sensitivity of imaging atmospheric Cherenkov tele- 
scopes decreases approximately with the square root of the 
solid angle of the gamma-ray emission region, and there- 
fore linearly with the source extension. For this reason the 

^ Quantities found by other authors using a different cos- 
mology are scaled to the value of .ffo adopted in this work 
throughout the paper 
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detectability of a source is proportional to its gamma-ray 
luminosity divided by its size . If it is assumed that 
the X-ray size and the X-ray brightness of a cluster to- 
gether are a measure of its gamma-ray flux, then Fx / Rx 
can be used as a figure of merit. It has to be noted that 
this selection procedure prefers galaxy clusters that host a 
so-called cooling core at their center. In a cooling core clus- 
ter, the central gas density is large enough that the radia- 
tive cooling time due to thermal X-ray emission is shorter 
than the Hubble time (see Peterson & Fabian 120061 for a 
review). This large density of target material is favorable 
for hadronic production of gamma rays. However, since 
only a small fraction of the total gas mass is contained 
in the cooling core, this will increase the total gamma-ray 
luminosity by only a modest amount. These selection cri- 
teria were applied to the galaxy clusters of the REFLEX 
survey (Bohringer et al. I2004p . which covers two thirds of 
the southern sky and contains 447 X-ray bright clusters. 
Based on this selection procedure, Abell 496 was found to 
be a promising candidate for H.E.S.S. observations. 

Abell 496 is a nearby (z — 0.033) relaxed cluster that 
features such a cooling core at its center. It also shows cold 
front substructures (Dupke & White I2003P . Cold fronts 
are sharp edges in the X-ray surface brightness that are 
inconsistent with the appearance of shock fronts. In the 
case of Abell 496, the cold fronts are likely to be caused by 
the oscillation of the central galaxy around an equilibrium 
position, since the chemical abundance of the ICM is flat 
across these structures (Dupke & White [2003p . The ab- 
sence of a jump in the abundance suggests that material 
at both sides of the front is of the same origin and hence 
the cold front is not caused by a sub-cluster merger. 

As a second target, a massive cluster was observed with 
a quite deep exposure. For selecting this target a differ- 
ent procedure was adopted than for Abell 496. Clusters 
were evaluated according to their accretion power, which 
scales with M^/'^ (M is the total mass of the cluster, see 
Gabici & Blasi [2D03l ISOOij) . Then the estimated accre- 
tion power was converted into an 'accretion flux' by de- 
grading it by the inverse square of the target distance. 
Finally the source extension was also taken into account. 
Since, as described above, the sensitivity of imaging at- 
mospheric Cherenkov telescopes decreases approximately 
linearly with the source extension, the 'accretion flux' was 
further scaled with the inverse of the cluster radius. The 
aforementioned criteria were applied to the X-ray selected 
galaxy cluster catalog of Reiprich & Bohringer (|2002p . 
Following this selection procedure the cluster Abell 85 was 
found to be a promising target. 

Abell 85 is a galaxy cluster with a complex morphol- 
ogy at a redshift of 0.055. In X-rays this object shows two 
sub-clusters merging with the main cluster (Kempner et 
al [20021 Durret et al. 12005]) . Additionally, Abell 85 fea- 
tures a cooling core at its center which is quite uncom- 
mon for merging clusters. Presumably the merging sub- 
clusters have not yet reached the central region of the 
cluster and have therefore not disrupted the existing cool- 
ing core (Kempner et al. I2002p . 



In the adopted cosmological model, the redshift of 
Abell 496 (Abell 85) corresponds to a distance of 134 Mpc 
(220 Mpc) and 1° relates to 2.35 Mpc (3.86 Mpc) at the 
target. 

3. Observations and data analysis 

The observations were performed with the H.E.S.S. tele- 
scope array, consisting of four lACTs located in Namibia 
(23° 16'18" S 16'^30 '00" E). The system is described by 
e.g. IHofmann 20031 Data were taken for Abell 496 from 
October to December 2005 and in October 2006. Abell 
85 was observed in October and November 2006 and in 
August 2007. The observations were performed in wobble 
mode (Aharonian et al. I2006P with the target being typ- 
ically at 0.7"^ offset from the center of the field-of-view. 
This allows a simultaneous background estimation in the 
same field-of-view. In total, 14.6 (32.5) hours of live time 
on Abell 496 (Abell 85) meet the standard data quality se- 
lection criteria and are used for the analysis. The average 
zenith angle of the observations of Abell 496 (Abell 85) 
was 28*^ (18°) resulting in a post-analysis energy thresh- 
old of 570 GeV (460 GeV) for hard cuts (see Aharonian et 
al. l2006|l . 

The data analysis after calibration consists of a recon- 
struction of properties of impinging gamma rays through 
the Hillas-parametrization of the observed Cherenkov im- 
ages (Hillas [T996|) . a background subtraction, and an esti- 
mation of the flux (Aharonian et al. I2006p . The direction 
of the primary gamma photon is estimated using geomet- 
rical analysis of multiple camera images according to the 
algorithm 1 as specifled by Hofmann et al. (|1999p . The 
energy of the primary gamma is estimated by comparing 
the image intensity with the values predicted by Monte 
Carlo simulations for the reconstructed gamma-ray direc- 
tion and impact point (see Aharonian et al. (|2006p for de- 
tails). A detailed comparison of the intensity of Cherenkov 
images generated by fast muons going through the tele- 
scope with simple theoretical calculation provides an ab- 
solute calibration of the optical efficiency of each telescope 
(Aharonian et al. l2006p . 

A system of cuts on parameters of each shower is used 
to reduce the prevailing hadronic background. Hard cuts 
from Aharonian et al. (20Q6P are used because they were 
optimized for sources of low fiux and for hard spectra. 

The remaining background is estimated using 
the reflected regions background estimation method 



(Berge et al. 20071. The method uses a number of control 
regions in the same field-of-view (FOV) as the observed 
target for estimating the background flux. In case of a very 
extended emitting region (radius > 0.7°) the background 
is estimated using the ON/OFF method (Weekes et al. 
I1989p . The background is then estimated from off-source 
runs. To assure an accurate background subtraction, these 
runs were performed under the same conditions as the on- 
source runs. 

All upper limits are derived using the approach of 
Feldman and Cousins (Feldman & Cousins I1998P at a 
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Table 1. Properties of the two galaxy clusters. Lx, 0.1-2.4 is the luminosity in the 0.1 - 2.4 keV band and R^^^ is the 
virial radius. Refs.: a) Reiprich & Bohringer (|2002p . b) Markevitch et al. (|1999p . c) Durret et al. (|2005p . 



cluster 


Mtot [Msun] 


Lx, 0.1-2.4 [erg s '] 


T [keV] 


Rvir [Mpc] 


redshift 


Distance [Mpc] 


Refs 


Abell 496 


3.1x10'* 


1.9x10*" 


4.7 


1.4 


0.033 


134 


a,b 


Abell 85 


7.6x10^'' 


4.8x10"* 


7 


1.9 


0.055 


220 


a,c 



99.9% confidence level. A power law gamma-ray spectrum 
with spectral index F = 2.1 is assumed. This value is based 
upon the expectations of Volk et al. p996p . In order to 
check the dependence of the results on the assumed index, 
the upper limits for F = 2.3 are also given. The difference 
is less than 10% in all cases. 

4. Results 

Abell 496 

Since the size of the VHE gamma-ray emitting region is 
unknown, four different analyses are used, integrating the 
signal over several sizes of on-source regions at the position 
of Abell 496. The integration radii 6 are chosen accord- 
ing to characteristic length scales of the density profile of 
the ICM, which acts as the target material for hadronic 
gamma-ray production. The density p of the ICM gener- 



ally follows a (3 model (Cavaliere & Fusco-Femiano 1976) 



Po 



1+ ( - 



21 -3/3/2 



(1) 



In the /3 model po is the central gas density, Tc is the 
core radius and the parameter (3 describes the slope of the 
density profile. 

The core analysis is used to search for a signal coming 
from the X-ray core region of the cluster. The core size of 
the P model was measured by Markevitch et al. (|1999p as 
178 kpc, corresponding to 0.08°. This size is comparable 
to the size of the point-spread-function of H.E.S.S. The 
signal is therefore integrated over a larger area with radius 
9i = 0.1°, optimized for a search of point-like sources 
(Aharonian et al. [2006)) . 

No significant emission is found from the central re- 
gion above the hard cuts threshold energy of 570 GeV. In 
order to search for other point sources off the center of 
Abell 496 a map was created with significances of point- 
like TeV gamma-ray signals on a grid around the center 
of the cluster (Fig. [1]). The significance distribution over 
the FOV is consistent with background fluctuations. The 
upper limit on the integral flux from the core region is 
derived to be F{> 570 GcV) < 4.8 x lO^^^ g-i ^^-2 
assuming F = 2.1 and F(> 570 GeV) < 5.2 x lO^^^ ph. 
s~^ cm~^ for F = 2.3, corresponding to 0.9 and 1.0% of 
the integral Crab Nebula flux above the same energy as 
measured by Aharonian et al. p006p . The integral upper 
limits as a function of the energy above which the flux is 
shown in Fig. [2] 



•2 -12''30 
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Fig. 1. Correlated significance map of the region around 
Abell 496 obtained using the core analysis. For each 
point in the map, the significance is evaluated by count- 
ing gamma-ray candidates within a radius of 0.1°. White 
dashed circles depict the on-source regions used for the 
core analysis, the IMpc analysis and the extended analy- 
sis (corresponding radii 6*1 = 0.1°, 6*4 = 0.4°, 6*2 = 0.6°). 
The very extended analysis is not shown since the on- 
source region is larger than the map. The distribution of 
significances in the FOV is consistent with background 
fluctuations. Also shown are black contours from hard- 
band ROSAT PSPC observations (Durret et al. 2000), 
smoothed by the H.E.S.S. angular resolution. 

The extended analysis is performed on the whole region 
of the X-ray emission from the cluster, using 02 = 0.6° as 
measured by Reiprich & Bohringer p002p . No signal is 
found above the energy of 570 GeV. Assuming a photon 
index of F — 2.1, an upper limit on the integral flux is de- 
termined F(> 570 GeV) < 2.4xl0~^^ph. cm^^, which 
corresponds to 4.6% of the integral Crab flux. Figure [2] 
compares this result with the core analysis. Assuming 
F = 2.3, the upper limit is F{> 570 GeV) < 2.6 x lO^^^p^. 
s-i cm-2 (5.0% Crab flux). 

The very extended analysis is aimed at investigating 
possible emission from the accretion shocks. A very large 
on-source region of radius 6*3 = 1.5° is used, corresponding 
to 3.5 Mpc. The data set is reduced to 9.8 live hours, be- 
cause not every on-source run had a corresponding good- 
quality off-source run that could be used for the back- 
ground estimation. Again, no significant signal is found. 
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Fig. 2. The H.E.S.S. 99.9% upper limits on the integral 
flux from Abell 496 as a function of the energy above which 
the flux is integrated, assuming four different sizes of the 
emission region. All curves are obtained by comparing the 
gamma-ray excess above a given reconstructed energy to 
the number of events expected for a spectral index of F = 
2.1. The result for the core analysis exhibits fluctuations 
caused by a lower number of excess events. Note that the 
upper limit for the very extended analysis was produced 
with a reduced dataset. 



An upper limit is determined to be F{> 570 GeV) < 
5.8 X 10-i2ph. s-i cm-2 (10.9% Crab flux) for T = 2.1 
and F(> 570 GeV) < 6.2 x lO^^^ph. s^^ cm^^ (11.7% 
Crab flux) for T = 2.3. 

An additional analysis was performed, using the inte- 
grating radius 64 = 0.4° that corresponds to a radius of 
1 Mpc from the center of Abell 496. The radius of 1 Mpc 
is not physically motivated, but other physical quantities 
that are used in the discussion are well measured within 
this radius. Hence these results are used for modeling in 
section[5] For the same reason, the upper limits are in this 
case calculated above 1 TeV. The resulting upper limits 
are F{> ITeV) < 7.2 x lO'^^ph. s'^ cm^^ for F = 2.1 
(3.2% Crab flux) and F(> ITeV) < 7.5 x lO^i^ph. 
cm-2 (3.3% Crab flux)for F = 2.3. 



All results are summarized in table [2l 



Abell 85 

Similarly to the case of Abell 496, several integration radii 
9 were probed to search for a signal from Abell 85. Abell 85 
could be observed at higher elevations which led to a lower 
energy threshold for this target. 

The size of the core region is 226 kpc (Mohr et al. 
I1999p . corresponding to '--^ 3'. The core analysis thus again 
uses 9i = 0.1°, the radius optimized for a search for point- 



like sources. No significant signal is found and an upper 
limit on the integral flux is derived F{> 460 GeV) < 3.9 x 
10~i3pj^^ g-i 55/^ (.^^-^^ fj^^-) Jqj. p ^ 2.1 and F(> 

460 GeV) < 4.1 x lO^^ph. cm^^ (q g% Crab flux) 
for F = 2.3. A significance map of the region around the 
center of the cluster is shown (Fig. [3|). Also in this case 
the significance distribution is consistent with background 
fluctuations. 
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Fig. 3. Correlated significance map of the region around 
Abell 85 obtained using the core analysis in the same way 
as in Fig. [1] White dashed circles depict the on-source 
regions used for the core analysis, 1 Mpc analysis and ex- 
tended analysis (corresponding radii 6'i — 0.1°, 6*4 = 0.26°, 
82 = 0.49°). The distribution of significances in the FOV is 
consistent with background fluctuations. Also shown are 
black contours from hard-band ROSAT PSPC observa- 
tions (Pislar et al. iT997|) . 

The extended analysis probes the region of X-ray over- 
density, as measured by Reiprich & Bohringer p002p to 
be O2 = 0.49°. The integral upper limit in this case is 
F{> 460 GeV) < 1.5 x lO'^^ph. s'^ cm^^ (2.0% Crab 
flux) for F = 2.1 and F(> 460 GeV) < 1.6 x IQ-^^ph. s^^ 
cm-2 (2.2% Crab flux) for F = 2.3. 

The very extended analysis uses an on-source region of 
radius 6*3 = 0.91° corresponding to 3.5 Mpc. The data set 
is reduced to 8.6 live hours due to a lack of appropriate off- 
source data. No significant signal is found. An upper limit 
is determined to be F{> 460 GeV) < 9.9 x IQ-^^ph. s^^ 
cm-2 (13.6% Crab flux) for F = 2.1 and F(> 460 GeV) < 
1.1 X 10-"ph. s"^ cm-2 (15.1% Crab flux) for F = 2.3. 

The next analysis probes the 1 Mpc region, corre- 
sponding to 64 — 0.26°. Again no signal is found and 
upper limits are derived: F{> ITeV) < 3.2 x lO^^^ph. 
s-i cm-2 (1.4% Crab flux) for F = 2.1 and F{> ITeV) < 
3.3 X 10-i3ph. s~^ cm-2 (1.4% Crab flux) for F = 2.3. 

The 95% X-ray containment (Perkins et al. I2006P re- 
gion is 510 kpc (based on parameters from Mohr et al. 
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Table 2. Summary of the results of Abell 496 and Abell 85 observations. 



Abell 496 



Analysis 


Radius 

n 


Radius 

[Mpc] 


Eth 
[TeV] 


Assumed 

r 


[10- 


Fui(> Eth) 
ph. cm"^ s--^] 


Fui(> Eth) 
[% Crab flux] 


Core 


0.1 


0.2 


0.57 


2.1 
2.3 




0.48 
0.52 


0.9 
1.0 


1 Mpc 


0.4 


1.0 


1.0 


2.1 
2.3 




0.72 
0.75 


3.2 
3.3 


Extended 


0.6 


1.4 


0.57 


2.1 
2.3 




2.4 
2.6 


4.6 
5.0 


Very extended 


1.5 


3.5 


0.57 


2.1 

2.3 




5.8 
6.2 


10.9 
11.7 


Abell 85 


Analysis 


Radius 

n 


Radius 

[Mpc] 


Eth 
[TeV] 


Assumed 

r 


[10- 


Fui(> Eth) 
12 ph. cm-2 s-i] 


Fui(> Eth) 
[% Crab flux] 


Core 


0.10 


0.4 


0.46 


2.1 

2.3 




0.39 
0.41 


0.5 
0.6 


95% X-ray containment 


0.13 


0.5 


0.46 


2.1 

2.3 




0.34 
0.36 


0.5 
0.5 


1 Mpc 


0.26 


1.0 


1.0 
1.0 


2.1 
2.3 




3.2 
3.3 


1.4 
1.4 


Extended 


0.49 


1.9 


0.46 


2.1 
2.3 




1.5 
1.6 


2.0 
2.2 


Very extended 


0.91 


3.5 


0.46 


2.1 
2.3 




9.9 
11.0 


13.6 
15.1 



I1999p . corresponding to an angular cut of ^5 = 0.13°. 
No signal is found and the upper limit on integral flux 
is F(> 460 GeV) < 3.4 x lO^i^pi^^ g-i ^^^^-2 g^/^ ^^.^^^ 

flux) for r = 2.1 and F{> 460 GeV) < 3.6 x lO-^^ph. s^^ 
cm-2 (0.5% Crab flux) for T = 2.3. 

The upper limits on the integral fluxes are plotted in 
Fig. m An independent calibration and analysis method 
combining a semi-analytical shower model and the Hillas 
analysis (Id e Naurois et al. 2005[) was used as a cross-check 
for both clusters, yielding consistent results. 

5. Discussion 

In the following, the production mechanisms of CRs are 
discussed and applied to estimate the total energy of non- 
thermal particles in the galaxy clusters (Section l5.ip . The 
upper limits in the TeV band obtained in the Section|3]are 
used to experimentally constrain the non-thermal content 
of Abefl 496 and AbeU 85 in Section O 

5.1. Estimates of energy in non-thermal particles 

For the understanding of the non-thermal history of 
galaxy clusters, the important parameter is the fraction 
of energy that is contained in non-thermal particles with 
respect to the thermal energy of the ICM (Eth). The to- 
tal energy of non-thermal particles is dominated by the 
hadronic component, since no losses are expected for par- 
ticles with energies of up to ~10^^ eV during the Hubble 
time. In contrary to that VHE electrons suffer from energy 



losses and do not accumulate in galaxy clusters. Therefore, 
for all the following calculations only hadronic CRs are 
considered. 

In clusters of galaxies CRs are accelerated by two dis- 
tinct mechanisms: externally via accretion shocks and hi- 
erarchical merger events, or via internal mechanisms such 
as supernova driven galactic winds and AGN outbursts. 
The energy in CRs produced by external processes will be 
proportional to the thermal energy of the ICM if it is as- 
sumed that the thermal plasma was shock-heated during 
the assembly of the cluster and that the CRs were acceler- 
ated in the same large scale shocks. The energetics of CRs 
resulting from supernova related processes will depend on 
the kinetic energy released by the supernova activity dur- 
ing the entire Hubble time, and thus on the total num- 
ber of supernovae exploding in the cluster volume (Volk 
et al. I1996P . Finally the contribution from AGNs will be 
constrained by the kinetic power of the outbursts and the 
timescale of activity. The total non-thermal energy is given 
by: 

-Enonth = ecxt-E-khi + ^int^'kin + -fAGN^active (2) 

In this equation, -Bnonth is the total energy in CRs, 
Eoxt is the efflciency of the accretion shock in accelerating 
CR particles, E^f^ is the kinetic energy of the accretion 
process which assembled the present day cluster, eint is 
the efficiency of supernovae and galactic winds in acceler- 
ating particles, is the combined kinetic energy pro- 
vided by all supernovae exploding in the cluster volume. 
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Fig. 4. The H.E.S.S. 99.9% upper limits on the integral 
flux assuming different sizes of the emission region. The 
curves are obtained in the same way as in Fig.O All curves 
were produced assuming a power-law with a spectral index 
of r = 2.1. The upper limit for ^5 = 0.13° is not shown 
for the sake of readability. It lies however very close to the 
curve of point source analysis. 

Pagn is the mean power of AGN activity and tactivo is the 
combined duration of all AGN outbursts. Values for the 
particle acceleration efficiency of shocks in galaxy clusters 
are uncertain and expected typical numbers are used for 
these estimates. The contributions to Snonth for the three 
described mechanisms are briefly estimated in the follow- 
ing. 

Accretion shocks: In galaxy clusters the accretion pro- 
cess is the dominant mechanism that heats the ICM. 
Intrinsic processes such as supernovae provide only a 
minor contribution to the thermal energy of the ICM 
(Renzini I2003P . Hence the thermal energy is a good mea- 
sure of the kinetic energy of the infall of the cluster 
building blocks. The thermal energy of the ICM can be de- 
termined by measuring the temperature and the density 
profile of the ICM. Here for Abell 496, a uniform tem- 
perature of 4.7 keV throughout the cluster volume and a 
[3 density profile for the ICM is adopted (numbers from 
Markevitch et al. |1999|) . A value of 4.3 x 10^^ gj.g ^j. ^j^g 
thermal energy is computed. For Abell 85, a temperature 
of 7 keV and the density profile of Pfrommer & Enfilin 
P004p is used. With these input values the thermal en- 
ergy of the ICM in Abell 85 is found to be 2.4 x 10*^^ 
erg. The value of eext should be close to 0.1 (e.g. Dorfi 
I199ip for strong shocks, but authors using nonlinear sim- 
ulations claim that it could be as large as 0.5 (Kang et al. 
I2002p . However, if the accreted material is already hot, the 
sound velocity in such a medium is large. Hence typical 
Mach numbers of shocks in such media are small and these 



weak shocks are less efficient in accelerating particles (e.g. 
Dorfi & V61k[T996|- This fact could reduce the value of ecxt 
in case of accretion of hot gas during sub-cluster mergers 
significantly. Keshet et al. (2004) have argued that only a 
fraction of 8 - 17% of the baryonic material of a galaxy 
cluster is accreted in strong shocks. This accretion is also 
the contribution to i?th that is accompanied by effective 
CR acceleration because the rest of the material is assem- 
bled in weak shocks with inefficient non-thermal particle 
production. Here it is assumed that 10% of the ICM is 
accreted in strong shocks with Ccxt = 0.3. Consequently, 
in this scenario the energy in the non-thermal component 
will be 3% of the thermal energy for both clusters. 

Supernova activity: To constrain the energetics of in- 
ternal processes in the cluster, the total cluster iron mass 
can be used to estimate the total number of supernovae 
(and thus the combined kinetic energy) that occurred in 
the cluster volume over a Hubble time. Supernovae come 
in general in two distinct classes which both feature about 
the same mean kinetic energy but differ considerably in 
the amount of iron produced: type la supernovae gen- 
erate on average 0.7 M© of iron whereas the iron yield 
of core collapse supernovae is 0.07 Mq (see Renzini 120031 
for a review). The ejecta of supernovae la can be distin- 
guished from the ejecta of core collapse supernovae by 
the ratio of iron to oxygen (Renzini I1993P . By knowing 
the relative contribution of supernovae la and core col- 
lapse supernovae, respectively, and using the mean iron 
mass produced by each class of events, it is possible to 
evaluate the total number of supernovae that have dis- 
tributed the observed iron. From the Fe/0 ratio found 
in Abell 496 and Abell 85 it is evident that type la su- 
pernovae contribute significantly to the chemical enrich- 
ment of the central cooling region (see Tamura et al. l2004p 
but are unimportant for the total iron mass in the entire 
cluster (De Grandi et al. I2004p . Hence it is assumed that 
the observed iron mass of 2.1 x 10^° M© (Abell 496) and 
5.0 X 101° Mq (Abell 85) (De Grandi et al. [2004)1 was 
entirely generated by core collapse supernovae, and there- 
fore an iron production mass of 0.07 M0 per supernovae 
(Renzini I1993P is adopted. As a result of these consider- 
ations it is found that 3 x 10^^ supernovae (Abell 496) 
and 7.1 x 10^^ supernovae (Abell 85) are necessary to en- 
rich the ICM to the observed level with iron. By adopting 
the canonical energy per supernova to be 10^^ erg, a total 
energy of internal processes of 3 x 10^^ erg (Abell 496) 
and 7.1 X 10^^ erg (Abell 85) is derived. This energy is in- 
jected into the ICM in the form of kinetic energy, thermal 
energy and CRs, and furthermore, depending on the envi- 
ronment of the supernova explosions, parts of it will also 
go to radiation losses (Dorfi , 1991, Thornton et al. I1998p . 
When assuming that 10% of the initial kinetic energy of 
the supernovae is converted into CRs (eint = 0.1 e.g. Dorfi 
[T99T|) either in the supernova shocks themselves or in the 
termination shocks of supernova-driven galactic winds, it 
is found that supernova activity produces a component of 
high energy particles with about 7% (Abell 496) and 3% 
(Abell 85) of the thermal energy of the ICM. Note that 
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eint can be smaller than 0.1 if a large number of super- 
novae explode directly in the hot ICM (Domainko et al. 
120041 Zaritsky et al. I2004p because in this case SN shocks 
are less efficient in accelerating CRs (Dorfi & Volk 1996]). 

AGN activity: AGNs are generally considered to be 
a major source of CRs in galaxy clusters (e.g. EnBlin et 
al. ri997 . AharonianHOni Hinton et al. [2UU7| . In coohng 
core clusters, the AGN activity is usually dominated by a 
powerful central AGN. In the case of Abell 496 no AGN 
- ICM interaction is observed (Dunn & Fabian 120061) but 
this could simply be due to a low activity period of the 
central galaxy at present. It has to be noted that in several 
galaxy clusters, buoyantly rising bubbles filled with radio 
emitting relativistic electrons are found as remnants of 
past AGN outbursts (e.g. McNamara et al. I200H Fabian 
et al. l2002p . which is also not seen in the cluster Abell 496 
(Dunn & Fabian [2006|) . The situation is different for the 
cluster Abell 85. In this cluster bubbles filled with non- 
thermal electrons injected by a past AGN outburst can be 
seen in the ICM. These bubbles are about 10^ years old 
and the energy that is necessary to infiate these bubbles is 
~ 10^^ erg (Dunn et al. 2005). In general, AGN can act in 
repetitive outbursts with up to 100 cycles over the lifetime 
of the galaxy cluster. In each cycle the AGN is active 
for ~10% of the time (McNamara et al. |200T1) . One may 
estimate the contribution of potential AGN activity in the 
past to the non-thermal particle component in this cluster. 
If it is assumed that a powerful AGN injects 10*^ erg/s 
of hadronic CRs in the ICM and is active for 10% of the 
cluster lifetime (which is essentially the Hubble time) then 
this mechanism will distribute about 3 x 10^^ ergs of CRs. 
Hence in an optimistic scenario, AGNs have the ability 
to provide non-thermal particles with an energy of 7% 
(Abell 496) and 1% (Abell 85) of the thermal energy of 
the ICM to the galaxy cluster. 

In summary, it is found that all three processes (ac- 
cretion shocks, supernova activity and AGN outbursts) 
may contribute a comparable share to i^nonth and it is 
further found that in total about 17% (Abell 496) and 7% 
(Abell 85) of the thermal energy of the cluster can be in 
the form of CRs. As stated above these numbers have large 
uncertainties and typical values are given. Additionally it 
is shown that the contribution to -Enonth by internal pro- 
cesses is more important for less massive, cooler clusters. 

5.2. Probing the non-thermal energy content of the 
cluster 

From the upper limits on the integral VHE gamma-ray 
flux obtained by the presented H.E.S.S. observations, it is 
possible to estimate the upper limit on the total energy 
of all hadronic CRs accelerated in these clusters within 
their lifetime for energies that are large enough for tt'^ 
production, but where particles are still confined in the 
system. 

The following estimates of the upper limits on the to- 
tal energy in hadronic CRs are given for a radius of 1 Mpc 



(0.4° for Abell 496 and 0.26° for AbeU 85) for both clusters 
(see Table[3]). For Abell 496, according to Markevitch et al. 
(fT999l) . the temperature of the ICM is 4.7 keV and the to- 
tal gas mass of the thermal ICM is 3.2 x 10"'^'^ Mq within a 
radius of 1 Mpc. The typical uncertainty of the quantities 
obtained from X-ray observations is about 10%. From the 
H.E.S.S. observations the upper limits presented in Sec. [4] 
are found in this volume. For Abell 85 a temperature of 7 
keV (Durret et al. [MIS]) and a gas mass of 1.2 x lO^'* Mg 
(following the density profile of Pfrommer & Enfilin [2004p 
is used. From the H.E.S.S. observations the upper hmits 
presented in Sec. |4] are found within a radius of 1 Mpc for 
AbeU 85. 

Since the distribution of CRs in galaxy clusters is not 
known, the implications of these upper limits on differ- 
ent scenarios of CR injection and propagation are investi- 
gated. First, a model with a constant CR density through- 
out the entire cluster is applied. This distribution is sup- 
ported by a scenario where CRs are accelerated in accre- 
tion shocks at the outskirts of the cluster and are mainly 
transported to the higher gas density regions of the cluster 
center by bulk motions of the ICM caused by hierarchical 
merger events (see e.g. Miniati 120031 for a discussion on 
the possible distributions of CRs in galaxy clusters). 

As a second scenario, it is assumed that the CR density 
follows the large scale density distribution of the thermal 
gas excluding the central cooling region. This is motivated 
by a scenario where the CRs are mainly injected by cluster 
galaxies which are more concentrated towards the cluster 
center. 

Finally, a somewhat more extreme model is adopted 
where the density in CRs in the central cooling region in- 
cludes an additional spike, similar to the gas density. It has 
to be noted that such a centrally peaked distribution of 
CRs should not form in clusters. Magneto-hydrodynamic 
instabilities are expected to lead to a de-mixing of gas 
and relativistic particles where relativistic particles with a 
softer equation of state will rise buoyantly to larger cluster 
radii and will show a less centrally concentrated distribu- 
tion than the thermal gas (e.g. Parker [T966l Breitschwerdt 
et al. 119931 Chandran & Dennis [2006p . Furthermore, since 
the central cooling region is much smaller than the whole 
galaxy cluster, CRs can leave this region also due to dif- 
fusion. 

In addition to the numbers computed for a radius of 1 
Mpc, upper limits of -Enonth/-Eth for the whole cluster are 
calculated with the upper limits on the gamma radiation 
found for a radius of the clusters that contains 95% of the 
expected gamma-ray emission (Perkins et al. 12006^ . The 
95% containment radius is computed with the assumption 
that the density in CRs follows the gas density including 
the high density cooling region (see Perkins et al. I2006p . 
The radius of 95% containment determined in this way 
is 1 Mpc for Abell 496 and 510 kpc for AbeU 85. Values 
of i?nonth/-Bth obtained with this specific method for both 
clusters are again given in Tab. [H It has to be noted that 
limits on -Enonth/£'th derived in this way are particularly 
small for clusters with a highly peaked gas density. This 
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results from deep upper limits on the gamma-ray bright- 
ness of the cluster due to a small spatial integration region, 
which can be adopted for the gamma-ray analysis for this 
case. However it has to be further noted that the value 
of £'noiith/£'th obtained in this way only holds for the as- 
sumption that the density of CRs is as peaked towards 
the cluster center as the gas density. This scenario seems 
to be disfavored by theoretical considerations as has been 
discussed in the previous paragraph. 

For the calculations of the upper limits of the non- 
thermal component of the cluster, the distribution of the 
thermal gas excluding the central cooling region (3' = 118 
kpc) found by Markevitch et al. (|1999p {(3 model) and 
the density profile obtained by Durret et al. (|2000p for the 
cluster Abell 496 including the cooling core (CC included) 
is used. In the case of Abell 85 the distribution of CRs 
is assumed to follow the single large scale f3 profile from 
Pfrommer & Enfilin (|2004p (/3 model) and to trace the 
combined double /3 profile with an additional component 
in the central high density cooling region (CC included) 
(again from Pfrommer & EnBlin [2004p . For computing the 
gamma-ray luminosities, the production rate of gamma 
rays per unit volume is integrated over the entire galaxy 
cluster assuming spherical symmetry. Calculations are ob- 
tained for the hadronic channel adopting the gamma-ray 
emissivity according to Drury et al. (|1994p . Results are 
given in Table [3] for a spectral index of the CR protons of 
2.1 and 2.3. A hard spectrum is expected in galaxy clusters 
since no losses of CRs should occur there, and therefore 
the observed spectrum should have the same energy dis- 
tribution as the primarily accelerated CRs. Depending on 
the model, hmits on E'nonth/^^th range from 0.40 to 1.12 
for AbeU 496 and from 0.03 to 0.15 for Abell 85 when a 
spectral index of 2.1 is assumed. 

Note that some models of other authors that concen- 
trate on external production mechanisms predict a ratio of 
CR energy to gas thermal energy of up to ~ 50% (Miniati 
et al. I200H Ryu et al. I2003p . Hence the upper limits con- 
strain models which favor a similarly large ratio of non- 
thermal to thermal energy. Limits obtained for Abell 85 
are especially interesting. These limits are for a hard spec- 
trum well within the prediction of the simple model devel- 
oped in Sec. l5.2l and exclude an unduly large component of 
hadronic CRs (i?nonth/£^th > 0.15). The non-detection of 
Abell 85 may even confirm the aforementioned theoretical 
arguments for which a very centrally concentrated distri- 
bution of CRs is disfavored. Indeed in this case, an upper 
limit of i?nonth/£^th < 0.03 can be derived which is smaller 
than model estimates (see Sec. 15. ip . However it is impor- 
tant to mention that for a steeper spectrum of the CR 
protons, the limit on the energetics of the non-thermal 
component would be larger than these optimistic model 
predictions. For stronger constraints on the component of 
non-thermal particles in galaxy clusters, longer exposures 
are required. In this context it should be noted that the 
sensitivity of H.E.S.S. scales with the square of the ob- 
servation time and therefore, in order to reach a twice as 



Table 3. Upper limits on the ratio of energy in the non- 
thermal component with respect to the thermal energy of 
the ICM (£^nonth/^'th) for different spatial distributions of 
the CRs in the cluster. All numbers are given for a radius 
of 1 Mpc with the exception of the numbers for the 95% 
containment radius (see main text). Particles of the non- 
thermal component have much larger energies (> 1 GeV) 
than can be obtained thermally and are responsible for 
the gamma ray production whereas the thermal compo- 
nent can be observed in X-rays. A value of i?nonth/-E'th > 1 
means that the energy of the non-thermal component ex- 
ceeds the energy of the thermal component. This config- 
uration is unrealistic since it would require that shocks 
are more efficient in accelerating particles than in heating 
up the shocked medium. Therefore, the presented obser- 
vations can not constrain any models for such a case. The 
uncertainties on the values given in this table is about 
40% due to the errors in the X-ray and gamma-ray obser- 
vations. 







Abell 496 


Abell 85 


spectral 


spatial 


-BnonthZ-E'th 




index F-, 


distribution of CRs 






2.1 


constant 


1.12 


0.15 




/3 model 


0.51 


0.08 




CC included 


0.40 


0.06 


2.3 


constant 


5.66 


0.75 




/3 model 


2.56 


0.40 




CC included 


2.03 


0.30 


2.1 


CC included 








95% containment 








radius 


0.40 


0.03 



sensitive limit on the energy in CRs as presented here, it 
is necessary to perform a four times longer observation. 

Very recently, also constraining limits of Enonth/Eth 
at the few percent level for the galaxy cluster Abell 521 
have been derived from radio observations (Brunetti et al. 
I2008p . With this approach the synchrotron emission of the 
secondary electrons produced by the proton - proton inter- 
actions is tested. The limits obtained in this way depend 
on the cluster magnetic field and are thus complementary 
to the limits derived from gamma-ray observations, which 
are independent of the magnetic field in the cluster. 

6. Conclusions 

No significant point-like or extended gamma-ray flux F(> 
570 GeV (Abell 496) and F(> 460 GeV) (Abell 85) has 
been found in the H.E.S.S. observations. Upper limits on 
the VHE gamma-ray flux from both clusters are derived. 
With the presented results, it is possible to constrain the 
energy fraction in hadronic CRs. For the case where the 
energy density of the CRs follows the density of the ther- 
mal gas and where the CRs have a spectral index of 2.1 
this is not more than 51% (AbeU 496) and 8% (AbeU 85) 
of the thermal energy. Especially for the case of Abell 85 
the upper limits obtained with H.E.S.S. can already con- 
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strain model predictions on i?nonth/-£'th- These limits are 
the best determined so far without a dependency on the 
magnetic field. However, for a uniform energy density of 
CRs or a softer spectrum of the particles, the limit on the 
fraction of energy in hadronic CRs versus thermal energy 
of the ICM is not constraining. Our theoretical consider- 
ations also suggest that the most promising clusters for 
VHE gamma-ray observations are those with the largest 
mass (Mtot > 10"'^^ ^q) and highest temperatures (T > 
7 keV) and, as far as supernova activity is concerned, 
those with the largest iron mass (Mpc > 5x10^" M0). 
Observations of this kind of galaxy cluster at a distance of 
~100 Mpc with the present generation of lACTs to a flux 
of F{> 1 TeV) < 10-12 pii^ g-i within a radius of 1° 

will provide the ability to test models of the non-thermal 
hadronic component in galaxy clusters, which predict a 
considerable fraction (^10%) of energy in CRs. 
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